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The G'2-MSSM is a model of particle physics coupled to raoduli fields with interesting 
phenomenology both for colliders and astrophysical experiments. In this paper we consider 
a more general model - whose moduli Kahler potential is a completely arbitrary G'2-holonomy 
Kahler potential and whose matter Kahler potential is also more general. We prove that the 
vacuum structure and spectrum of BSM particles is largely unchanged in this much more 
general class of theories. In particular, gaugino masses are still suppressed relative to the 
gravitino mass and moduli masses. We also consider the effects of higher order corrections 
to the matter Kahler potential and find a connection between the nature of the LSP and 
flavor effects. 
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I. INTRODUCTION 

From several theoretical points of view, the existence of moduli fields seems inevitable. For instance, 
supersymmetry may be the mechanism responsible for stabilizing the scale of the Standard Model. Su- 
persymmetry requires supergravity, whose only (known) reasonable UV completion seems to be String 
theory; and along with string theory come extra dimensions and their moduli. In fact, since string/M 
theory contains no dimensionless parameters, moduli appear necessary to explain the observed values of 
various couplings in nature. From the bottom up, moduli appear in various theories with "dynamical 
couplings" as well as in Inflation - the inflaton field is usually a neutral scalar field aka a modulus. For 
all of these reasons and more, moduli physics and phenomena must be considered seriously. 

In a series of papers, 01, 0]) 0], a very detailed model of moduli physics coupled to matter has 
been described. The G2-MSSM model, largely inspired by M theory compactifications on manifolds 
of G2-holonomy, is a model in which strong gauge dynamics in the hidden sector generates a potential 
which both stabilizes all the moduli fields and simultaneously generates a hierarchically small scale - thus 
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solving (most of) the hierarchy problem. The model has an interesting spectrum: moduli have masses in 
the 50-100 TeV region, scalar superpartners and higgsinos have masses in the lO's of TeV region, whilst 
gauginos, which are the lightest BSM particles have masses of order lOO's of GeV. Direct production of 
gluinos and electroweak gauginos are the dominant new physics channels at the LHC. The nature of the 
LSP is also very interesting as it is a neutral Wino. Moreover, its production in the early universe is 
dominated by the decays of the moduli fields (ie non-thermal production) and can naturally account for 
the observed fraction of dark matter today. The moduli and gravitino problems are avoided due to the 
gravitino mass scale being one to two orders of magnitude larger than the TeV scale. One drawback of 
the model is the fine tuning between the lO's of TeV scale and Mz and is the reason the model solves 
most of the hierarchy problem and not all of it. 

However, the G2-MSSM model, as defined in Q], is based on some specific assumptions about the 
moduli and matter Kahler potentials, albeit with the claim that these are general enough to incorporate 
all of the essential ingredients of more general Kahler potentials (and hence G2-manifolds) . Thus far, 
there has been no serious study of these assumptions and it is the main aim of this paper to undertake 
this. The main result that we prove here is that the mass spectrum of the theory depends very weakly 
on the specific form of the moduli Kahler potential; in fact the spectrum depends on the Kahler potential 
for moduli only through the fact that it is the Log of a homogeneous function (the volume of the extra 
dimensions); the precise nature of this homogeneous function is fairly irrelevant as we will see. 

We also discuss the Kahler potential for charged matter fields. We give three consistent arguments 
for calculating the moduli dependence of the matter kinetic terms in 4d Einstein frame. Whilst non- 
trivial, these modifications do not change the results of [H, [3] much. More importantly, we also consider 
higher order terms in the matter Kahler potential, in particular the terms which are usually considered 
troublesome for flavor physics in theories of gravity mediated susy breaking. Whilst we expect that such 
operators will be suppressed, if they enter with large coefficients they can affect the mass spectrum: 1) 
they can directly alter the scalar and higgsino masses, which are typically large; 2) they can indirectly 
(via threshold effects from the higgsinos) alter the nature of the LSP. In particular, we find that the LSP 
can also be a Bino in some cases. This provides a connection between flavor physics and the nature of 
the LSP in models of this sort. 

The paper is organized as follows. The next section describes some simple properties of the moduli 
space metric for general G2-manifolds which will be important for our later considerations. Section III is 
devoted to the Kahler potential for charged matter fields. Following this, we re-do the analysis of Moduli 
stabilization from [1] in this much more general context. In section V we compute the mass spectrum and 
susy breaking couplings in the minimum of the potential and demonstrate that it is almost identical to 
that of the original G2-MSSM. In section VI we present a further generalization of the construction. In 
section VII we renormalize the Lagrangian down to the Electroweak scale and give the spectrum there. 

II. GENERAL PROPERTIES OF MODULI SPACE METRICS ON G2 HOLONOMY 

MANIFOLDS 

In this section we describe some very general and simple properties of the moduli space metric of 
G2-niaiiifolds. It is these simple properties, which will allow us to draw very general conclusions. 

The metric g{X) on a G2 holonomy manifold X can be expressed in terms of the associative three-form 
$ as 



gij = (dets) 9 Sij 



(1) 



with 




klnmrst pl2...7 



+ 1. 



(2) 
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Expanding $ in terms of basis harmonic three- forms (pi G H'^{X, Z) (modulo torsion) we obtain 



$ = ^ Si(j)i , N = h^{X) = dim {M{X)) , (3) 

where Sj are geometric moduh corresponding to the perturbations of the internal metric. The complexified 
moduli space M {X) of a G2 holonomy compactification manifold X has holomorphic coordinates given 
by 

Zi = ti + isi , (4) 

where ti are the axions parameterizing the zero modes of the 11-dimensional supergravity three form C3. 
The classical moduli space metric (not including possible quantum corrections) can be derived from the 
following Kahler potential [^] 

i^ = -31n47r^/Vx, (5) 

where the dimensionless volume Vx = Vol{X)/l\.^ is a homogeneous function of Sj of degree 7/3 and Im 
is the lid Planck length. The homogeneity of Vx is the key property that we will utilize in what follows. 
In terms of the associative three- form the volume is given by [J] 

Vx = ^ y $ A . (6) 

Define the following derivatives with respect to the moduli 

dk , d^k 
—— and Kij = 

OSi OSiOSj 



ki = — and ki, EE ^j-^ . (7) 



The matrix kij , the Hessian of k is related to the actual Kahler metric Gjj which controls the kinetic 
terms as 4(5 = , where in the Hessian we simply replace index j with j . Since Vx is a homogeneous 
function of degree 7/3, the first derivative of k defined above has the following property 

N 



Y^Siki = -7. (8) 



i=l 

Differentiating ^ with respect to sj we obtain an important property of the metric Kij 

N N 

Sikij = —kj , and since kij is symmetric ^jkij = —ki . (9) 
i=i j=i 

Let us now introduce a set of dual coordinates {rj} defined by 

dVx 



dsi 



(10) 



Note that the variables {rj} are homogeneous functions of {si} of degree 4/3 . Using the homogeneity of 
the volume together with the definition (fTO|) we can express the volume Vx as 

7 ^ dVx ^ 3 ^ 

-Vx = Y,Si-Q^ = Y,SiTi, y^ = -^s.ri, (11) 

i=l 1=1 1=1 
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and 



dK 3 dVx 3Ti , , 

dsi Vx dsi Vx 



Combining ([3]), ([6]) and ([IT]) we can reexpress the dual variables as 



Ti = l [ (t>i/\*'^ = l [ , (13) 

where for each harmonic basis three- form (pi G H^{X,Zi) we introduced a Poincare dual four-cycle [tj] S 
H4^(X). We now use the above duality to make a particularly convenient choice of the basis harmonic 
three- forms. In particular, we choose a basis {4>i} S H^{X) such that the periods of the fundamental 
co-associative four-form over the Poincare dual four cycles are positive definite: 

*$>0. (14) 
This choice of a basis becomes obvious when we recall that for a generic basis four-cycle [tj] G H/i{X) 

V0l{[Ti]) >\f *^ = Ti, (15) 

where the above relation becomes an equality if and only if the corresponding four-cycle is co-associative. 

All geometric moduli describing the fluctuations of the internal metric must be massive in order to 
satisfy constraints from fifth force experiments and cosmology. At the same time the vacuum expectation 
values of the moduli (coordinates Si on the moduli space) must be fixed in the region of the moduli space 
where the geometric description makes sense. In Section IIVI we describe a way to stabilize the moduli, 
which ensures that we find isolated minima that satisfy these conditions automatically. 

It turns out that for our purposes it is convenient to introduce a set of "angular" variables a-i defined 

by 

Oj = SiKi = -4— , no sum over i . (16) 

3 Vx 



We see that a, are scale-independent and satisfy 

N 



Ea. = |. (17) 



i=l 



Thus, we can also parametrize the moduli space M[X) by a subset of A^ — 1 variables ai plus one volume, 
e.g. the volume of the manifold Vx- Differentiating the Oj allows us to introduce the matrix 



which has components 



dai 

Pij = —Sj-z — , no sum over j . (18) 

OSj 



Pij = ^^ijSjKi + SiSj^kij , no sum over i, j . (19) 



Pij has the following contraction properties, which follow from (jl7p and the fact that Oj are homogeneous 
of degree zero 

N N 

^Pi,=0, and ^Pi,=0. (20) 
i=i j=i 
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We can then write 



where the matrix Aij is defined as 



K,, = ^A,j , (21) 



Ai, = + ^ , (22) 



and satisfies the following contraction properties 

N N 

Ajj = 1 , and ^ AijGj = Ui , (23) 
i=i j=i 

where we used (j20p to derive (|23p . Note that parameters Oj defined in (jl6p are the components of an 
eigenvector a of the non-Hermitian matrix A with unit eigenvalue. 

We can compute the formal inverse of the Hessian metric, K^^ . By definition of the inverse it must 
satisfy 

N 

Y,K''K,k = Si, (24) 



i^^i = £i£i(A-i)^i, (25) 



and using (|2T|) it can be expressed as 



where the inverse matrix (A satisfies 



Symbolically we can express A ^ as 



N 

Y,{^-'f^jk = Sl. (26) 



= — p , (27) 

a 



which in terms of components translates into 

{A-^yi = 6'^ - P^j - + Pa-Pij - - Pii-Pim—Pmj -+.... (28) 

ttj ai Qj ai Qm O-j 

Using (|2Up and (|28p we derive the following properties of the inverse matrix 

AT N 

^{A~^f = 1 , and ^{A~^y^aj = at , (29) 
i=i i=i 



which could have also been obtained directly from (j23p . Note that although we do not have a closed 
form expression for the components (A~^)*-^, the contraction properties in (j29p are what will ultimately 
allow us to derive explicit expressions for the terms in the soft breaking lagrangian - since such couplings 
depend only on the contractions and not the precise details of the functional form of Vx- Before going 
on to the details of these calculations, we first must consider the Kahler potential for matter fields in M 
theory. 
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III. KAHLER POTENTIAL FOR CHARGED CHIRAL MATTER 

In this section we re-visit the Kahler potential for charged matter fields in M theory. In practice, the 
absence of a useful microscopic formulation makes it difficult to compute the moduli dependence of the 
Kahler potential for these fields in general. Below we outline three arguments for the structure of the 
Kahler potential - first from dimensional reduction, second based on the scaling properties of physical 
Yukawa couplings and the third based on the form of the threshold corrections to the physical gauge 
coupling. Happily, all three methods agree. 



A. Kahler potential from dimensional reduction 

In M theory, charged chiral matter is localized near conical singularities 0], 0, 0], @]- These are 
literally points in the seven extra dimensions. Because of this, we expect that the kinetic terms for the 
chiral matter fields should be "largely independent of bulk moduli fields" that the G2 manifold X has. 
The precise meaning of this statement will be clarified below in terms of the scaling property of the kinetic 
term. They could, of course depend on local moduli inherent to the conical singularity, but, since, in a 
supersymmetric theory, a single chiral multiplet in a complex representation of the gauge group usually 
has no D or F-fiat directions [9], there are typically no such local moduli. 

There is a subtlety in the above general arguments. Since, in four dimensions, a scalar field kinetic 
term is not invariant under Weyl rescalings of the metric, one has to pick a Weyl gauge. We will argue 
that the correct Weyl gauge for the statement above is NOT the 4d Einstein frame. Therefore, the kinetic 
term for chiral matter will be non-trivial in the 4d Einstein frame, which is the standard one in which to 
define the Kahler potential. 

Since the physics of a conical singularity in M theory does not introduce any new scale, asides from 
the lid Planck scale, the only reasonable Weyl frame is the lid Einstein frame. Therefore the lagrangian 
density in the lid frame is 

L ~ Mn^^i? + 5^^ dM(t)dN4>'^ 9^''' K{si)^i + . . . , (30) 

where 5j is a delta function peaked at the position of the matter multiplet containing the scalar field (j) 
and has mass dimension seven. n{si) is a homogeneous function of the moduli of degree zero which will 
generally be of order one and vary adiabaticallyHJ, i.e. 

fc=i 

The above property implies that K{si) remains invariant when the moduli are rescaled as Si — >• Asj, 
thus explicitly implementing the idea that in the lid frame the kinetic term of a matter field localized 
at a point p ^ X \s "largely independent of bulk moduli". A particularly simple example satisfying ([3T|) 
is when n{si) = const. Integrating this over X leads to a 4d density 

L4 ~ VxMl^^^Ri + Kisi)g^''d^(f>d,(l>^^4. , (32) 

where Vx is the volume of the extra dimensions in lid units. This is the Lagrangian in lid Einstein 
frame. If we now Weyl rescale into the 4d Einstein frame we find 

L4 ~ ^V9eRe + '^^^g^Ed^,(t)du(l)^^/g^ , (33) 



^ In the toroidal Type IIA compactifLcations with intersectiong D6-branes, the analog of k is a homogeneous function of 
the complex structure moduli of degree zero since it only depends on the relative angles Of [loj . 
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where the subscript E indicates that we are using the 4d Einstein frame metric. 

We have only considered the Einstein-Hilbert and kinetic terms of the matter fields. Including all the 
other terms would give the 4d supergravity Lagrangian in Einstein frame. In particular, from this we 
would read off that the Kahler metric for the multiplet containing (j) is 

If we introduce dimensionless fields (p as in (p = mpicj), the Kahler potential is 



K = K{si)^ = Hi{si)^mli . (35) 
As we will see, this is consistent with the arguments given in the next subsection. 



B. Kahler metric from the properties of the physical Yukawa couplings 
Here we will describe an alternative way of deducing the volume dependence of the Kahler metric 



for charged chiral matter. This method is due to Conlon, Cremades and Quevedo [ill] and utilizes 
the relation between the physical (normalized) Yukawa couplings l^/s-y and the unnormalized Yukawa 
couplings Y^p^ that appear in the supergravity superpotential. Recall that in G2 compactifications of 
M theory, a superpotential Yukawa coupling Y^^^^ between the multiplets a, /3, 7 that are localized at 
three co-dimension seven singularities is induced by an M2-brane instanton wrapping a supersymmetric 
three-cycle connecting the three singular points. The absolute value of the Yukawa coupling is given by 

^ e-''-'''^^-' , (36) 

where 

Vap-y = ^ ml'^^Si , (37) 

i 

is the volume of the supersymmeric three-cycle. After we diagonalize the Kahler metric for the mat- 
ter fields and go to the canonical basis, the relation between the absolute values of the physical and 

unnormalized Yukawa couplings is simply a rescaling by e^^'^ (^KaKpK^ 

\Yo.p^\ = e"^/' \YU^ iko^KpkX^'^ ~ \Y'\ (v^K^kpkX^'^ . (38) 



On the other hand, one can construct perfectly well-defined seven-dimensional local models where 
the G2 manifold is non-compact, e.g. an ALE-fibration over a three-sphere or a quotient thereof, in 
which case Vx — >• 00 =^ TUpi/Mn — )• 00 and gravity is effectively decoupled. Such models can also contain 
charged chiral matter fields and since their interactions are determined locally, the corresponding physical 
Yukawa couplings should not vanish when gravity is decoupled. Local models of this type can be obtained 
via lifting effective theories on intersecting D6-branes in Type IIA to M theory. 

Therefore, locality implies that the physical Yukawa couplings should be independent of the overall 
volume Vx in the limit Vx — >• 00. For that to happen, the Kahler metrics K^, Kp, in (|38p must scale 
with the volume Vx as 

k^r^kp^k^^ X^ (39) 



which is in perfect agreement with the form of the Kahler metric derived in the previous subsection. 
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C. Consistency check for the Kahler metric 



In this section we confirm the form of the Kahler metric for charged chiral matter by comparing 
the threshold corrections to the physical gauge couplings in G2 compactifications of M theory with the 
general results in M = 1 D = A super gravity. 

Let us first consider a hidden sector containing a pure glue SU{N) supersymmetric Yang-Mills theory. 
Using the notation in we have the following relation for the gauge coupling at one loop 

167r2 167r2 / ^ 

where S are the one-loop threshold corrections, g{fi) is the physical gauge coupling and gM is the tree-level 
Wilsonian gauge coupling. In our convention qm is related to the gauge kinetic function / as 

^=Im/. (41) 
9 m 

Recall that the Wilsonian gauge coupling gets renormalized at one loop only. On the other hand the 
physical coupling g(fi) is renormalized to all orders. The threshold corrections come from massive states 
and are independent of the scale //. Based on the topological arguments [l^l, the threshold corrections 
due to Kaluza Klein modes in G2 compactifications of M theory are rather simple and can be calculated 
even without knowing the G2 metric! Friedmann and Witten explicitly computed one-loop threshold 
corrections due to the heavy Kaluza-Klein modes living on a supersymmetric cycle Q with 61 (Q) = 
^2(2) = and a non-trivial fundamental group. Such corrections come in a form of linear combinations 
of Ray-Singer analytic torsions [1^] . which are topological invariants of Q. For the case at hand, the 
threshold corrections are given by 

S = 2Nln VqA^ + 2 ^ TiTin, , (42) 

i 

where Vq is the volume of the supersymetric cycle Q, 71 are the Ray-Singer torsions corresponding 
to different irreducible representations of the fundamental group and Q are the generators of SU{N). 
Here, the cutoff dependence appears as a correction due to the zero mode contributions transforming in 
the trivial representation of vri(Q). Once the threshold corrections are included explicitly, a somewhat 
unexpected cancellation of the A-dependence occurs pjj] and the one-loop relation can be written as 

l"'^ l"'^ ,mn(']^S- (43) 



where 

S' = 2^TlTrn.Q\ (44) 

i 

In fact, the cancellation of the A-dependence occurs for any supersymmetric cycle Q with 61 (Q) = 

62(2) = 0. 

Now we would like to consider a more general case when the gauge theory is a supersymmetric 
QCD with Nf flavors of chiral matter fields Qa transforming in of SU{N) plus Nf flavors of Qa 
transforming in N. Each chiral matter field transforming in a complex representation arises from a 
separate co-dimension seven conical singularity on X with each singular point Pi G Q. It was argued in 
[12l | that the singularities producing charged chiral matter fields have no effect on the KK harmonics of 
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the seven-dimensional vector multiplet. Moreover, since the conical singularities introduce no new scale 
below the eleven-dimensional Planck scale Mn, the effective cutoff scale for these multiplets is naturally 
Mil. Including such multiplets into the running is straightforward and results in 

In addition to the KK thresholds, there may be some unknown corrections due to possible charged 
massive matter fields with masses of order Mn. At this point we cannot say with certainty whether 
such massive charged M theory modes are present in the spectrum but we cannot exclude this possibility 
either. Just like the KK thresholds, these corrections cannot be holomorphic functions of the chiral 
multiplets Zi describing the moduli of X since the axion partners of the geometric moduli decouple from 
the computations of the threshold corrections. However, there may be some non-holomorphic as well 
as constant contributions from such massive charged states. For now we will simply assume that they 
are constant and result in a slight shift of the tree-level gauge coupling. On the other hand, moduli 
dependent contributions may arise from non-perturbative corrections due to membrane instantons but 
they will be exponentially suppressed and can be safely neglected. 

Our next task is to independently verify that the Kahler metric for the charged chiral matter fields 



matches the previously obtained result (j34p . Here we will use a strategy similar to the one in 1J|, 1C| 
and compare ()45p with the corresponding one- loop expression in M = 1 D = 4 supergravity given by 15 1 

167r^ 167r^ / "i", \ . / I 



(37V - Nf)ln - {N - Nf) K + 2N In (^-J-^ - 2Nf In (^K^a) ■ (46) 



In the above expression, K = —SIuAtt^^^Vx is the Kahler potential for the moduli and Kaa is the Kahler 
metric for the charged chiral matter fields. We can use the definition of the four-dimensional Newton's 
constant K4 = ^/SttGn = l/nipi in terms of the eleven-dimensional gravitational coupling kh 



K 
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Vxll 



(47) 



M 

in combination with the common convention 2k^^ = (2tt)^M^^ and Mn = 271 /Im to obtain 

Ml = . (48) 

Using the above relations together with At: / g\j = Vq/1\^ we have from ([l6l) 

.2 1 R^2 



= ^ - (3/^ - iV/) In ( - 2iV, m 4.V3 + 2iV In ( j^) - 2N , 1„ (a'^ 



= ^-3]Vln^ ' j +iV,l„(^) -2A';l„(vi/f„,) -l„((2,r(8.r0 . (49) 

The appearance of the last term is most likely due to the convention used to define Mn in terms of kh 
as well as the ambiguity in defining the relation between Im and Mn. Thus, we shall regard this term 
as an artifact and ignore it in further discussion. Comparing (j49p with the expression on the right hand 
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side in (j45p we conclude that up to a constant multiplicative factor, Kahler metric for the charged chiral 
matter fields Qa is 

~ , (50) 

vx 

which precisely matches the result obtained in the previous subsections. On the other hand, the constant 
term S' in (j45p has no corresponding analog in (j49p and represents a genuine threshold correction to the 
Wilsonian gauge coupling qm- 

In the framework of = 1 D = 4 supergravity, the RG-invariant scale where super QCD with > Nf 
becomes strongly coupled is 

where the second exponential factor is due to the local SUSY. The Affleck-Dine-Seiberg effective super- 
potential [3] W should be identified with 

■iN-Nf 

(52) 

det(QQ)'^-^/ 

where the gauge coupling inside A is complexified. Using (|5ip . up to an overall numerical constant we 
obtain 

W ^{N- Nf)m^^^'''f det(QQ)"^^e'^^e"^(^^ . (53) 
In (j53p . the dimensionful chiral matter field^ Q can be expressed in terms of dimensionless fields Q as 

Q = 'mpiQ , and Q = rripiQ . (54) 

Then, the superpotential becomes 

W = C{N - Nf) mlidei{QQ)'^ e~^^^^ , (55) 

where C is an overall numerical constant. In our further notation, we also define the following constants 

s' 

C = C e~^^^^ , and A = {N-Nf)C. (56) 
Let us now consider the case of Nf = 1 flavors. Introducing an effective meson degree of freedom 



'2QQ, (57) 

we can rewrite the superpotential in terms of (j) as 

W = Am^i(j)~^e^^ , (58) 

where we have absorbed the factor of 2^/^^"^^ into the normalization constant C. Along the D-flat 
direction we have Q = Q and the Kahler potential for the matter fields can be rewritten in terms of the 
effective meson fields (f) as 

K = n[.Si)-^mpi + K{si)-^mpi = ^(si)— m^, . (59) 



^ Here we suppressed the flavor index. 
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D. Higher order terms 



Based on three independent arguments we have been able to deduce the volume dependence of the 
Kahler metric for charged chiral matter fields localized at co-dimension seven singularities. Denoting the 
visible sector charged chiral matter fields by Q'^ their Kahler potential is then given by 

;f = M(f!)«!5l!. (60) 

Vx 

In the regime where the size of the super symmetric cycle supporting the visible sector is large (this 
assumption is justified in the context of the MSSM where the corresponding volume is chquj' ~ 25 ) we 
can perform a systematic expansion in the inverse volume of the cycle (weak coupling) so that in the 
leading order K,^p{si) is a homogeneous function of Si of degree A, satisfying 

^i—Q^. = ^f^a^iSi) ■ (61) 

i=l * 

Based on the property that a given charged chiral matter multiplet is localized at a point p £ X we 
expect that A = 0. Nevertheless, for the sake of generality we will keep A as a free parameter for the time 
being. In our derivation of the Kahler potential we so far neglected possible higher order contributions 
to the visible sector matter Kahler potential of the forrcH 

6K = c^p{s,)Q:QfM, + ... = c^p{s,)91^^ + . . . , (62) 

which gravitationally couple the hidden sector meson to the visible sector fields Qa. In the above 
expression, the subscript c denotes canonically normalized matter fields in the 4-dimensional Einstein 
frame. Such couplings can create problems if the meson i^'-term is quite large (which is true in the G2- 
MSSM) because they can induce flavor changing neutral currents. This is the flavorproblem of gravity 
mediated susy breaking mo delfl These terms were neglected in our previous work [2[]. 

Technically, computing the unknown coefficients c^p{si) from the underlying theory is difficult, goes 
well beyond the scope of this work and our aim here is not to explain the flavor structure of the super- 
symmetry breaking Lagrangian. Rather, we would like to understand the effect that the presence of such 
terms might have on other sectors of the theory, e.g. their effect on superpartner masses and couplings. 
For these purposes it is sufficient to assume that the flavor structure of the Kahler metric is completely 
determined by the matrix K^p[si), so that 

c(£i) 
3 

where we introduced the factor of 1/3 for future convenience. As we will see in the later sections, whilst 
this does not introduce any flavor violation, the point will be that the effect of such terms on the mass 
spectrum will be similar even if we introduced flavor violating terms, as should become clear even tuallyS 
Actually, such a form might arise from an expansion of the Kahler potential if the visible and hidden 
sectors were completely sequestered. Though we do not expect M theory to be sequestered, it can be 
useful to think of the sequestering as an extreme limit in a more general model. 



c^-p{si) = n^-p{si) , (63) 



Here we set rUpi = 1 and treat all matter fields as dimensionless in units of mpi. 

Note that the Kahler metric derived from (jGOp introduces no flavor problems since, as we will see from explicit com- 
putations, the mass matrix for the scalars will be proportional to Hafti^i) and therefore diagonalization and canonical 
normalization of the kinetic terms automatically results in universal scalar masses. 

Generically, the absence of flavor changing neutral currents implies that the off-diagonal entries in the mass matrix for 
the canonically normalized squarks and sleptons are suppressed, though in particular models even stronger constraints 
are necessary, e.g. the requirement that the diagonal entries are nearly degenerate, depending on the spectrum and the 
A-terms [Hi. 
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A sequestered Kahler potential has the form 



K--? = -3 In ( - - -At„^Q"Q^^ ) , (64) 



—d>d) / 

3^^ 3 

and the Kahler metric for the visible sector is given by 

47rV3y^ _ i<^<^ ■ ^^^^ 

Absorbing the factor of Att^^^ into the definition of the fields and expanding the above expression in 
powers of (p'^/Vx we obtain 

K7 = ^fl + 4^') + ---- (66) 



Comparing the above expression with (|62p we can read off the coefficients 



CI = l^'c.B , (67) 



which corresponds to (|63p when c(si) = 1. Hence, function c{si) in (|63p is the measure of deviation of the 
matter Kahler potential from the exactly sequestered form. As was pointed out in [l^, sequestering is not 
at all generic in string/M theory and presumably G2 compactifications of M theory are no exception. 
We thus will regard the value of c{si) in a given vacuum as a parameter and consider the theory for 
various values of c(sj). 

Combining all of the previous considerations, the visible sector matter Kahler metric and its inverse 
take the form 



Ks = ^ l + cis.)^ , (68) 



k^'^ = K^\s,)-, ^^-T^ « n'^'\s,)Vx[ 1 - c{s,) 



where K°'^{si) satisfies 

K''^{s.i)Kp^{si) = 6^. (69) 

Combining (jOip with the above we conclude that K°'^{si) is a homogeneous function of the moduli of 
degree —A. Function c{si) will be typically assumed to take values in the range 

< c{si) < 1 . (70) 

However, as long as the Kahler metric is positive-definite, one may also consider the regime when c{si) < 0. 
Diagonalizing the Kahler metric of the visible sector we obtain 

Kad^P=i^UKypU^p, (71) 

where the eigenvalues are given by 

i^a = ^fl + c(..)#V (72) 



Vx \ ^ "^Vx 
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and Haisi) are homogeneous functions of degree A that satisfy (j6ip . In computing the anomaly mediated 
contribution to the gaugino masses, it wih be necessary to compute various derivatives of InKa- For this 
purpose, it turns out that it is very convenient to express Ini^'o, as 

InKa = In Ka{si) - In Vx + In ^1 + c{si)-^-^ » + hiKa(si) + (c(si) - 1) + const , (73) 

where K is the Kahler potential in 



IV. MODULI STABILIZATION 

In this section we reconsider the problem of moduli stabilization with the much more general moduli 
and matter Kahler potentials introduced in the previous section. We will be working in the frame- 
work of = 1 D = 4 effective supergravity and will demonstrate that all the moduli can be stabilized 
self-consistently in the regime where the supergravity approximation is valid. Recall that in the compact- 
ifications we study here, non-Abelian gauge fields arise from co-dimension four singularities [l^, 21] 



[23 |. [H]- In other words, there exist three-dimensional submanifolds Q inside the G2-manifold X, along 
which there is an orbifold singularity of A-D-E type. The basic idea is that strong dynamics in the hidden 
sector breaks supersymmetry, stabilizes the moduli and generates a small scale. The simplest possibility 
consistent with the supergravity approximation is a hidden sector with two gauge groups SU{P + Nj) 
and SU{Q) where the first is super QCD with Nf = 1 flavor of quarks Q and Q transforming in a 
complex (conjugate) representation of SU{P + 1) (the corresponding associative cycle Q contains two 
isolated singularities of co-dimension seven) and the second hidden sector with the gauge group SU{Q) 
is a "pure glue" super Yang-Mills theory. One can easily consider more general gauge groups without 
much qualitative difference. One can also consider a setup with charged matter in both hidden sectors. 
However, as was demonstrated in [ij, in such cases, one of the two F-terms coming from the matter 
fields in the hidden sectors is always suppressed relative to the other and thus does not contribute to the 
quantities relevant for phenomenology. A single hidden sector gauge theory is also enough to stabilise 
the moduli, though the vacuum is not in a place where supergravity is trustable! 

Therefore, the non-perturbative effective superpotential generated by the strong gauge dynamics in 
the hidden sectors is given by 

W = Ai(j)^e''''^ + A2e'''^^ . (74) 



The matter field (p represents an effective meson degree of freedom defined in (I57p in terms of the chiral 
matter fields Q and Q. The coefficients 6i, 62 and a are 

27r ^ 27r 2 ,^ , 

h = -, h = -, a = --. (75) 

In [ij it was explained that if one uses a superpotential of the form (j74p . de Sitter vacua arise only 
when Q > P (if we include matter in both hidden sectors dS vacua exist without such condition). Hence, 
we will keep this in mind from now on. 

In ([7^) we explicitly assumed that the associative cycles supporting both hidden sectors are in pro- 
portional homology classes which results in the gauge kinetic function being given by essentially the same 
integer combination of the moduli for both hidden sectors 

N 



f = Y^NiZi, (76) 



i=l 
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were 

N 



lm{f) = VQ^ [ ^ = J2N^Si (77) 

is the volume of the corresponding associative cycle with the integers Ni specifying the homology class. 
While one can certainly consider possibilities where the gauge kinetic functions /i and /2 are not pro- 
portional, the results in [l[ taught us that unless /i oc /2 it is more difficult to stabilize all the moduli in 
the regime where the supergravity approximation is valid. Thus, obtaining solutions which we can trust 
is the main reason for choosing to consider the case where /i = /2 = /• Obviously, progress in the more 
general cases would be welcome. 

Typical examples for 3-cycles supporting non-Abelian gauge fields in G2-manifolds are spheres and 
their quotients such as Lens spaces considered in 12]. The expression in (|74p can in principle 

contain many additional non-perturbative contributions if X contains other rigid associative cycles. In 
that respect, the two terms included in ([71}) should be regarded as the leading order exponentials. As 
long as Q and P are large enough compared to the Casimirs from the other gauge groups, the remaining 
terms will be exponentially suppressed in general. This is particularly true for the membrane instanton 
corrections to (I74p which come with exponentials containing bi = 2-k. On the other hand, some such 
instantons induce Yukawa interactions among the visible sector matter fields and are therefore implicitly 
assumed to be part of the full superpotential. 

The total Kahler potential - moduli plus hidden sector matter, is given by 

K = -^\nAT:^/'^Vx + ti{si)^ . (78) 

vx 

In what follows we first consider a simplified case where the function K{si) is a pure constant, i.e. 

K{si) = 1 . (79) 

However, in section IVll we will generalize our results to the case where ^(sj) is a homogeneous function 
satisfying (131 p . The important point is that even then the functional form of the soft breaking terms 
remains virtually unchanged compared to the simplified case, thus validating our approach. In general, 
()78p must include the contributions to the Kahler potential from all matter sectors including the visible 
sector as described in the previous section. However, since the visible sector fields will obtain zero vacuum 
expectation values (vevs), they can be dropped for the purposes of stabilizing moduli. 
The standard M = 1 D = A supergravity scalar potential is given by 

y = [K^'^FnFrn - 3|Ty|2) , (80) 

where the F-terms are 

Fi = diW + WdiK = iNie'^''^-^-hiAi4>le-^^^-' + b2A2e-^''^A (81) 



2si V ^Vx 

F^ = d^W + Wd^K = -e'''^^-^^~'^aAi(l)^"^e~'''^-^ (82) 



(83) 



Vx 

In the above we used 

dK 1 dK 



dzj 2i dsi 
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together with the definition of in ()16p in combination with 

d 1 Ki 



(84) 



We also parametrized the meson field <j) as 

<P = <Poe'' , (85) 

and fixed one combination of the axions and the meson phase 9 

cos{{bi- b2)N -t + aO) = -1. (86) 

Before we proceed to constructing de Sitter vacua it is instructive to take a step back and consider 
a simpler case where the first non-perturbative term in the superpotential is also a pure gaugino con- 
densate arising from a "pure glue" supersymmetric Yang-Mills theory. In this case one possible solution 
corresponds to a supersymmetric AdS extremum described by the following set of equations 



which is equivalent to 

3ai 



Si 



(-biAie-b^^-^+b2A2e-^2^-^ 



Using the contraction property (jl7p we can find from ([88]) that the volume • s of the hidden sector 
three-cycle can be determined by solving the following transcendental equation 



7 ( -Aie"''i^-^V A2e"''2^-' 

N-s = — -j^ (89) 

2 (-biAie'>'^^-^+b2A2e~''^^-^ 



In the limit when N ■ 1, the approximate solution is given by 
Vq = N-s^ , ln(44^)>0, when bi > b2 Sz Aibi > A2b2 , or bi < b2 Sz Aibi < A2b2 . (90) 

bl — 02 V^202/ 



The moduli vevs can then be found from 



where the seven-dimensional volume is stabilized at 

7/3 



yx = v'/'(^^ Vxisk) . (92) 



at. 



Note from (19ip that at the extremum, the periods of the co-associative four-form Tj ~ Ni up to a positive 
constant. Recall that here, Ni are the integers representing the homology class of Q: 



Ni= ^i, <PieH'{X, Z), (93) 
'[S] 
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where the harmonic three-form (pi is Poincare dual to a four-cycle [rj]. Therefore, from (j9ip we see that 
extremization of the supergravity scalar potential dynamically stabilizes the co-associative four-form 
to be proportional to the integral homology class of the associative three-cycle Q: 

Fi = =^ *<l> = a-PDxiQ), 0<a£R. (94) 

Therefore, in the basis specified by (|14|) the integers Ni must be positive definite 

A^i > 0, \fi = l, ... N . (95) 

In order to determine the values of Oj at the minimum we substitute our expressions for Si ()9ip into the 
definition of Oj in (|16p to get a system of N transcendental equations, which then completely determine 
Oj in principle 



+ 3Ni = . (96) 



Note that the dependence on Vq in pUj) is gone due to the scaling property of the volume Vx- Hence, we 
have recast the problem of determining the moduli vevs at the minimum into a problem of determining 
the values of aj. Obviously, obtaining general analytic solutions for from (I96p is impossible in practice, 
since Vx has not been specified. However, precisely because the moduli vevs at the minimum are given 
by (|9ip . it turns out that in order to compute the quantities relevant for particle physics, one does not 
need to know the values of Oj explicitly. All one actually needs to know are the contraction properties 
(dZ]) and (1291). 

Therefore, the results we derive will be valid for any singular manifold of G2 holonomy containing an 
associative three-cycle Q that contributes to the non-perturbative superpotential in a form of at least two 
gaugino condensates, whose integral homology class in the basis |j^[ ) is specified by positive integers. By 
explicitly checking in explicit toy examples, both numerically and analytically it seems that, for a given 
form of Vx, an isolated solution indeed exists. 

In principle, there exists an alternative way of determinimg Oj more directly, although in the long run 
it may be more practical to solve the system ()96J) . Namely, suppose one can reexpress the volume Vx 
in terms of the dual variables defined in (jlOp rI With respect to Tj the volume Vx is a homogeneous 
function of degree 7/4. Then, we find 

dVx _ ^ dVx dsj _ ^ dVx dsj _ dsj_ _ 3 

Qt- - 1^ f)^^ f,-^. - 2^ f)^. /)^. - Z^'^^f)^. - zl**' ^ > 



^ dsj dTi ^ dsj dTj ^ c?r,- 4 



where we used the property that Sj are homogeneous functions of Tj of degree 3/4 and the symmetry of 
the Jacobi matrix 

dr^ ^ d^Vx ^ 5r^_ ,gg. 
dsj dsjdsi dsi ' 

Then, using (|16p . ([U7|) and (pi]) we obtain 

s^n 4r, dVx 4^ 9 In Fx 



Vx 3Vx dn 3 dr. 



(99) 



One needs to ensure that in the new basis the signature of the Hessian ? Y,^ remains Lorentzian. 
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where in the final step we used the property that the "angular" variables Ui do not scale. Then we can 
re-express the moduli vevs ([9T|) as 



4,, dlnVx 



7 ^ dr, 



(100) 

T^=Ni 



In order to illustrate how the system ()96p is realized in practice we give a couple of explicit examples, 
though we stress that we have checked many more general examples than just those given here. Let us 
first consider a particularly simple A^-parameter family of Kahler potentials consistent with G2 holonomy 
where the volume Vx is given by 

N N 

yx = J{ s^l^ , where ^ = 3 • (101) 
1=1 1=1 

In this case the solutions to (|96p are simply constants given by 

ai = ni. (102) 

In fact, this example represents the class of Kahler potentials considered in the previous work [l] and the 
solutions are discussed in detail there. 

One may consider more complicated examples such as 

^ „fc ^7 
Vx = Vk , where Vk = Ck TT s"' , such that VA: nf = - . (103) 

o 

k i=l 1=1 

In this case system ()96p translates into 

T.{-'--^)-^Il[^) =0- (104) 

k j=i ^ 



In these cases one can check numerically that, for very generic sets of parameters {n^, c^, A'j}, the sys- 
tem of equations (I104p yields positive solutions for Oj, where the Hessian matrix H(Vx)ij has Lorentzian 
signature 3. For example, choosing A^i = 1, A'^s = 1, -^3 = 1, -^4 = 1 we numerically compute for the 
following toy examples with four moduli 

Vx = sfsjsfsf - -sjs|44 - -sf444 ^ 1.038 , 02 ~ 0.648 , 03 ^ 0.324 , 04 « 0.324 , 

3 2 

Vx = sfs^sf + \s\slsl + \sfslsl ^ ai^ 0.051 , 02 ~ 1.478 , 03 ^ 0.459 , 03 « 0.344 , (105) 



where for both examples 

d'Vx 



sign 



dsidsj 



(+,-,-,-), (106) 



which explicitly demonstrates that having positive solutions for is fairly generic and more importantly 
is guaranteed when Vx is not just a randomly picked homogeneous function of degree 7/3 but represents 
an actual volume of a G2 manifold X. 



This is necessary if the homogeneous function Vx is the volume of a genuine G2-manifold. This also guarantees positive 
kinetic terms for the moduli fields. 
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We now go on to consider de Sitter vacua by including the charged chiral matter fields Q and Q into 
the hidden sector. The superpotential and the Kahler potential are given by (j74p and (j78p . In order to 
compute the scalar potential we need to compute the inverse Kahler metric. Using the Kahler potential 
(|78p together with (jl6p . (j2ip . ()83p and ()84p we first obtain the following components for the Kahler metric 

Note that on the right hand side of the above expressions aj and Ajj are the same real quantities defined 
previously with index j replaced by j . 

The inverse Kahler metric must satisfy the following set of equations 

K^^ K-^^K^~^K-^^ = ^^, (108) 

After a little bit of work we obtain the following components for the inverse Kahler metric 

, ^ 4^ (A-lJ^ , ^ 2 ^ , _ -i^^ , (109) 

3Vx 3Vx SVx 

Note that despite the fact that the matter part of the Kahler potential in ()78p is only given up to the 

quadratic order in we decided to keep all the higher order terms inside the inverse Kahler metric. 

.2 

This is self-consistent as long as the combination appearing in the inverse Kahler metric is stabilized 
at a value sufficiently smaller than one such that the quartic and higher order terms are suppressed. 
Now, putting all the pieces together we obtain the scalar potential 



V 



647rV3 



|EEf^ii^^(M.*---^-M.e---f (110) 

. «=1 J =1 ^ ^ 3Vx 




3i I '/'o 3Vy 
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To understand the minima of the potential we will use the techniques developed earlier in Namely, 
we will work in the regime when the volume of the hidden sector associative cycle Vq = N ■ sis large and 
expand our solutions in the inverse powers of this volume. This is equivalent to an expansion in the UV 
weak hidden sector gauge coupling. In this long a tedious procedure we utilize the methods developed in 
[l|, yet with some important modifications. 

Since we are considering the simplified case by setting n^Si) = 1 in the Kahler potential for the effective 
meson field, the supersymmetry breaking F-term contributions are functions of Vx and N ■ s only and 
therefore the scale invariant "angular" coordinates ai will remain the same as in the supersymmetric case. 
On the other hand, the "radial" coordinate parametrized by Vq (or Vx) will be shifted. Reintroducing 
the notation of [ij 

^^A^^-(bi-b2)N-s^ y^bia-b2, z^bja-bl, (111) 

we therefore make the following ansatz for the moduli vevs at the minimum 

s^ = ^-L. (112) 

In this notation, the volume of the associative cycles supporting the hidden sector gauge groups is given 

by 

VQ = N-s = -LY,a, = --L, (113) 
in which case the moduli ansatz (I112p can be rewritten as 

■n = ^lv,. (114) 

Let us first assume that L is non-zero and finite when y ^ 0. This assumption will be verified in this 
section by determining L explicitly. Then, we get from ()113p and the definitions above 



Vq 



oo 



y^O 



a 



02 Vq 



(115) 



Vq = N-s 



■In 



62^2 



1 PQ 
2ttQ-P 



In 



QA 



i9o 



PAr, 



This fixes the value of the volume Vq of the hidden sector three-cycle. 

We now go on to demonstrate that the ansatz for the moduli vevs (I112p indeed represents the correct 
solution at the minimum of the scalar potential. In particular, we must verify our assumption that L is 
non-zero and finite in the limit y — )• by determining L self-consistently in this limit. Hence, we will 
now derive the equation for L and demonstrate explicitly that one of the possible solutions is indeed 
non-zero and finite in this limit. After minimizing the potential with respect to the moduli Sj and using 
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the definitions piip we obtain the following system of equations 



dV 



4 5 
3 5sfc 



1 + 



(116) 



3Vx 



3Vx 



1 + 



3l + 



3Vx 



^ ^ ^Vx 



3Vx Sk 



Ofc 7 



1+ 



3Vx 



3Fx SfcA^fe 3 



aa + 77— X 



dk Vx 



+— 1 + - 



7 1 



-2iVfc 



3 1 + ^ 3yx 

3Vx 



k2 \2 



7afc 



1 



9 gfc I 0^ 

3Vx 



a2 \2 



aa + 77— X 



1 + 



7 1 



3 1 + ^ 3yx 

^ 3Vx - 



aa + 77— X ) I aoia + — —x + r-—y 



Vx 



Vx SkNk Vx • 



A2e 



0, 



where in one of the intermediate steps we simplified 



Vx 



1 + 



^Vx 



3VxJ 



\ 2 



7 1 



<P ak 



Sk 



31+^ 3Vx Sk 
^ 3Vx 



1 + 



3Vx 



(117) 



Multiplying (jll6p by and using the explicit expression for the potential (jllOp in terms of the 
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quantities (jllip we obtain 



3 1 + 



- y ^ + 4 • 5 - + 7 1 + 



3Vx 



(118) 



aa 



2 \ 2 



+ 



3 ak dsk 



_J_l/ 



3 



a,: 



3Vx «J 



1 



x2i + 4L afe 



(-\ , -/-^ ^'3Fx 



.NkSky ,NkSk 



-(.N ■ s)^- - 4^(7V • s)- - 7-^ 



+2 



A^fcSfc / 2 
Ofc 1 3 



• s) '—yr- + -{N ■ s) 

3i + ^ 3^ l + JL 3^ ^ 



y/x (j)'^ ak 7y 



\ 



3Vx 



1 + 



3Vx 



3Vx SkNk 3 X 



aa (p 



(N-s) 



y/x '7 I / biaa cfP' y 



l + ^^3 V 



Vx 7 1 



,2 



3 1 + ^ 3Vx 

^ ^ 3Vx . 



X ~'~ Vx X 



+ 



V 



V ^ 3Vx 



_2^1x^^7 1 
afc (/>2 I 



oa / gfcia 0^ Qfc ^^'1 _ n 

3i + _|L3yx / V a: V ^ VxSkNk Vx x ^ ~ 



At first sight it appears that finding an analytic expression for L from (jllSp is hopeless since a closed form 
for (A"^)*-' is unknown and Oj have not been determined explicitly. However, upon further examination 
we notice that in order to find L from pisp we only need to know the contraction rules p7p . (j20p and 
([22D. Indeed, using the the ansatz (|112p together with the definition ([T5D and applying (fT7|) . (pJ|) and 
(l29l) we first evaluate the terms 



ik 



yZ 



(119) 



(120) 



and then use the same ansatz pi2p and contraction identities for the rest of the terms in pisp to obtain 
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the following equation for L 

^2 



3 1 + 



14 L 



1 + 



+ 1^ + 7(1 + 



7 



9 1 + ^ + 3 



+ 



56 zx -L'^ 



aa (p^ \ Vx 



31+^ 9y2-,,^ 



+ 



28 



^2 
L2 



1 + 



^Vx) 
7 1 



3l + ^3F.^ 



1 + 



3Vx 



3l0f 



-4L 



28 



28 



Wx 



~^Vx, 



:i2i) 



-2L 



+2L 



14 L 



3 1+^ 9 1 _L ^ y2 
^ ^ 3Vx ^ SVx 



14 

+ 



9 + 



3yx 



14 L 



1 + 



h — 

y Vx 



7 1 



/ aa 
' — + 



a2 \ 2 



3 1 I ^ 3yx / V 2: Vx 

^ 3Vx / ^ 



-2L 



16f 



1 + 



7 1 



aa (p 



-2 1 + 1 



3Vx 
a2 



3 1 + ^ 3Vx 

3Vx 



1 + 



3Vx 



abia (f) 



y 



+ 



Vx 



aa dp' , 
X Vx , 



aa 
X Vx 



aa 

— + 



a2 \ 2 



Vx 



Multiplying the above equation by ^(1 + ^^^) ™ and taking the limit y — )• we obtain 



aa\ 

3^y 



, _ ^2 ^iaay _ ^biaay 



3xz 



2xz 



1 + 



^ 3biaay f-^.W 
3Vx J Uxz V 9 Fx 



where we dropped terms of 0{y^) and higher. A non-trivial solution can be obtained by solving the 
corresponding quadratic equation 



3x 



^haay _ biaay ^ ^ _^ 



3xz 



2xz 



3Vx 



+ 



3biaay 
lixz 



1 + 12 ii 

9 Vx 



which is analogous to the equation in the second line in (126) of [l| 
Solving (I122p to the first subleading order in y results in 



aaVx 
(l>oX 



+ 1 



0(122) 



,^ , aa\ 3biaa 



1 I aaVx 
^ + —TIZ- 



Uxz 1 - H 



1 + 



3Vx 



(123) 



Hence, we see that this solution is non-zero and finite when y — )• and therefore is self-consistent. This 
is the solution describing the minimum of the potential. We must note that there is another possible 
solution of (jl22p for which L ~ y — t- 0. In fact this other solution corresponds to the extremum at the 
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top of the potential barrier and we will not discuss it further. Using (jl23p we can now compute the first 
subleading order correction to a to obtain 

P 7P(3(Q-P)-2) J_ 

12vrQ Vq ^^^^^ 

P 7{Q- P)2 2 \ P 



Q ' 2Q2 3{Q-P)J Peff ' 



where we have introduced 

P„„ = Pin I 1 

Using (|124p we can express the solution for L from (|123p as 



W^Pmrn^l. (12.) 



+ 1+ .„^''i „ H 1 + 



2ft// V (Q-P)*o/\ 3Vx 
In the leading order, the moduli vevs are given by 

_ tti SQPeff 

~ WiU7T{Q-P) ■ 



(127) 



We note that since Oj, Ni are positive, we need Peff > if Q > P, so that there exists a local minimum 
with Si > 0. 

The next step is to determine the vev of the effective meson field by minimizing the potential with 
respect to cpQ. Let us first compute the potential at the minimum as a function of the meson. The result 
is given in equation ()140p and the reader not interested in its derivation may proceed directly there. 
It turns out that since the moduli vevs at the minimum are proportional to Ui/Ni as in pi4p . explicit 
computation of the P-terms at the minimum and various contractions thereof while using the rules (jlZp 
and ()29p becomes possible. Let us demonstrate some of these computations in detail. First we need to 
identify the gravitino mass in terms of our notation in (jllip . Using the usual definition in combination 
with (jllip we have 

^3/2 = e^'/^lWl = e''/^\x\A2e-'^^-\ (128) 

Because the existence of de Sitter vacua requires g - P > (see m for details) we obtain using (|124p 
that 

x^^-KO. (129) 

On the other hand, since 7713/2 > we can express the following combination in terms of the gravitino 
mass 

e^/'^xA2e-^^^-' = -my2 ■ (130) 
We now multiply Pj in (jST]) by e^/^ and using (jllip and ()130p express 
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where denotes the overah phase of the superpotential. Using the ansatz ()112p for Sj we obtain from 

dm]) 



3Vq\ ^ 2 V 3Vx 



^3/2 
"^3/2 



where in the second hne we used 



''l = Ivq, (133) 



y 7 

obtained from (jll3p . Similarly, we find from (|8ip using (jllip together with ()130p 



^/2F, = e^(---^)f^ + #V3/2. (134) 



Before computing e^/'^F^ we would like to express the K^^ components of the inverse Kahler metric at 
the minimum using the ansatz (I112p for sj as follows 

3a, y 3a,N. ^ , J±_ "^^la^N- i , ^ ' ^ ' 

Contracting (jl32p and p34p with the inverse Kahler metric and using the solution for L from (I123p we 
then obtain 

^K/2pi ^ f,K/2j^ij p ^ e^/^K'^Fr (136) 



-te 



-i-yw 



■ 3 



1 _|_ Q^yx 

1 - If "^^/^ ^ 7v7l 

where in the last line we used (|124p to plug into x, y, and z defined by (jllip except for the combination 

^1 + "^a^^ ^ and kept the leading term in 1/Pe//- Note that in order to get from the second to third line 

in (1167p we used the second contraction property in (j29[) . Similarly, contracting (jl32[) and ()134p with the 
corresponding components of the inverse Kahler metric ()109p we obtain 



e 



Using the results (fT32]l . (fTSill . (fT36|) and (fT371l together with (fT23]l and (fTMD we can compute the 
following contributions 
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where we also used N ■ s = Vq while performing the computations in the first line of (|138p . 
Then, the potential at the minimum is given by 



Vo = e^{F'Fi + F^F^ - 3|Ty|2) (139) 



X 



eff 



+ ^\ ^rni,^ + —^{— + ^] (l + -^)(^) m2,^-3m2 



Using (jl24p and dropping the terms of order 0{1/ P"^^^) we obtain the following expression for the leading 
contribution to the vacuum energy as a function of the meson field 



Q-P VxJ PeffK 3{Q-P)J\Q-P VxJ Vx 



^^3/2 • (140) 



The polynomial in the square brackets in (jl40p is quadratic with respect to the canonically normalized 
meson vev squared (j)'^ = (^"q/Vx with the coefficient of the {4)q/Vx)'^ monomial being positive (+1) and 
therefore, the minimum Vq is positive when the corresponding discriminant is negative. Tuning the 
cosmological constant to zero is then equivalent to setting the discriminant of the above polynomial to 
zero, which boils down to a simple condition 

14(3(0 -P)- 2) , 
" 3(3(<3 - F) - 2V6(<? - P)) 

Note that Peff defined in (|125p is actually dependent on (p but because of the smallness of a and the Log 
dependence, it was safe to use the approximation Peff ^ const. This approximation turned out to be 
self-consistent since Peff is fairly large. From (I14ip we see immediately that 

P^jf >0 ^ Q-P>3. (142) 

Minimizing (jl40p with respect to (j)^ and imposing the condition that the expression in the square brackets 
in (jl40p is tuned to zero, we obtain the meson vev at the minimum in the leading order 

If we tune the leading contribution to the vacuum energy and set Q — -P = 3 we obtain 

^e//«63.5, :^~ 0.75. (144) 



Recalling the factor of two in the definition of the meson field (|57p we note that along the D-flat direction, 
the bilinears of the canonically normalized charged matter fields that appear in the original Kahler 
potential have a somewhat smaller vev 

< QcQl >=< QcQl >=< QcQc >~ 0.37 m^; , (145) 

which makes it a bit easier to justify the truncation of the higher order terms in the Kahler potential for 
hidden sector matter. 

We find numerically that for the minimum value Q — P = 3, the tuning of the cosmological constant 
by varying the constants Ai and A2 inside the superpotential results in fixing the value of Peff at 

61.648, (146) 
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while the canonicahy normahzed meson vev squared is stabilized at 

</,2 = ~ 0.746, (147) 
Vx 

thus confirming the analytical results above. For example, we obtain the values in (jl46p and (jl47p by 
minimizing the scalar potential numerically for the following toy example with two moduli 

P = 27, Q = 30, ^1 = 27, ^2 = 2.1544, Ni = N2 = I ,Vx = sfsl + \s\sl + \slsl 

^ si^ 34.52, S2 ~ 63.13. 

It is instructive to compare the moduli vevs obtained above numerically with the values obtained by 
using the analytic expression (jl27p . However, before we can apply (jl27p we need to determine the values 
of Oj at the minimum. This can be done by plugging the values of Ni, and Ck into the system (I104p 
and solving it numerically. To compute Si from (|127p we use Pg// from (jl46p in order to get a better 
accuracy. As a result, we obtain: 

aiss 0.825, 02 -1-51, siss34.7, §2 ~ 63.4, 

which confirms explicitly that the analytic expression (jl27p for the moduli vevs at the minimum is indeed 
very accurate and reliable. Here we also verified that the Hessian of the volume has Lorentzian signature. 

Although the value in ()147p is not much smaller than one, the combination inside the inverse 
Kahler metric (jl09p has a value 

'^^^ 0.25, (148) 



^Vx 

which is small enough to make the quartic and higher order terms which we kept inside the inverse Kahler 
metric much smaller. 

As we will see in the computations that follow, the value of Peff will enter into many quantities 
relevant for particle physics, such as tree-level gaugino masses, etc. Here we note that small changes 
in Peff do not lead affect the supersymmetry breaking masses much, but do change the cosmological 
constant significantly. For instance, while changing the value of Peff in the range 61 < Peff ^ 62 hardly 
affects the values of the soft breaking terms, as will be evident from the corresponding explicit formulas, 
such small changes in Peff result in vastly different values of the vacuum energy: 

61<Pe//<62 ^ -(m3/2mpi)'x 10-3 <yo< + (m3/2mp;)'x 10-3. (149) 

Therefore, once we coarsely tune Peff to its approximate value, the cosmological constant problem 
becomes completely decoupled from the rest of particle physics. Even though this should be the case, it 
is satisfying to see it explicitly in a complete example of moduli coupled to matter. 

Note also that in the original paper [l| we obtained Peff ~ 83. This is due to the different matter 
Kahler potential considered there. As we will see, this numerical difference will result in slightly different 
values for the soft breaking terms if compared to those obtained in [l|,[3]- 

Recall that for a stable minimum to exist it is necessary that Q — P > 3. We have seen that when 
Q — P = 3 and the minimum of the potential is approximately tuned to zero, the value of Peff ~ 60 which 
ensures that the moduli (jl27p can be reliably fixed at values large enough to satisfy the supergravity 
approximation. On the other hand, when Q — P = 4, from (I14ip we get Peff ~ 20, in which case if Q is 
fixed, the moduli vevs become smaller by about a factor of four. Thus, unless the ranks of hidden sector 
gauge groups are incredibly large, situations when Q — P > 3 may put our solutions well outside of the 
supergravity approximation. Therefore, from now on we will only consider the case when Q — P = 3 to 
ensure the validity of the regime where our construction is reliable. 



28 



V. MASSES AND SOFT SUPERSYMMETRY BREAKING TERMS 

A. Gravitino mass 

In supergravity the bare gravitino mass is defined as 

"^3/2 = mpie^/^\W\ = e^/Vl^2e"^2VQ ^ (^5q) 
and can now be computed since we stabilized Vq explicitly. It is given by 

When the cosmological constant is tuned such that ()146p is satisfied, for Q — P = 3 we obtain 

m3/2«9x 105(TeV)-^, (152) 

where C2 = A2/Q was defined in (|56p . Calculating C2 goes beyond the scope of this paper. Here we will 
treat C2 as a phenomenological parameter with values C2 ~ 0(0.1 — 1) since it may experience a mild 
exponential suppression as in ([56j) . 

On the other hand, the actual value at which the volume Vx must be stabilized can be almost uniquely 
determined from the scale of Grand Unification. In particular, we can use equation (4.12) in [l^ to express 

" 8^ " 327r2M2^^yx ' ^^^^^ 
where the factor L{Q) is due to the threshold corrections from the Kaluza-Klein modes and is given by 

L(Q) = 4gsin2 {5irw/q) , (154) 
such that 5w is not divisible by q. For typical values 

aGUT = ^ = ^, Mgj/t = 2 X 10^6 GeV , (155) 
Vq 

we obtain 

Vx = 137 A X L(Q)2/3 . (156) 

In Table I we list a few typical benchmark values for the volume and the resulting gravitino mass up 
to the overall factor C2. 

Interestingly, the gravitino mass scale naturally turns out to be constrained to 7713^2 ~ C'(10)TeV. 
While this is presumably large enough to alleviate the gravitino problem, it is also small enough to give 
some of the superpartners masses which can be easily accessible at the LHC energies. As we will see 
below, this is possible because of the significant suppression of the tree-level gaugino masses relative to 

"l3/2- 

In addition to the gravitino mass it is instructive to compute the scale of gaugino condensation. Using 
(jl46p the volume of the hidden sector associative cycle for Q — P = 3 is given by 

From (|5ip the scale of gaugino condensation in the second hidden sector is 

e-H^a £-20/3 11 X 10"GeV , , 
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Point 1 


Point 2 


Point 3 


Point 4 


Point 5 


Point 6 


Point 7 


q 


2 


3 


4 


4 


6 


6 


6 


w 


1 


1,2 


1,3 


2 


1,5 


2,4 


3 


Vx 


549.6 


594.5 


549.6 


872.4 


453.7 


943.7 


1143.2 


"^3/2/(^2 


70TeV 


62TeV 


70TeV 


35TeV 


93.0 TcV 


31 TcV 


23 TcV 



TABLE I. Typical values of Vx and 771,3/2 divided by C2 for different values of q and w. 



B. Moduli masses 



In order to compute the masses of the moduli we first need to evaluate the matrix Vmn with m, n 



1, + 1, with components given by 

^ O O ' ViN+l = r, r,, , Vn+1N+1 = r,, r. , , (159) 

OSiOSj OSiOcpo O(poocpo 

at the minimum of the potential. However, because the Kahler metric in (jlOTp is not diagonal, we 
also need to find a unitary transformation U which diagonalizes the Kahler metric. We denote all the 
components of the Kahler metric as Kmn- Then, by diagonalizing Kmn we obtain 

Kk5M = ul^K^nU^j. (160) 

After that, we need to rescale the fluctuations of the moduli around the minimum by the corresponding 
l/V 2Kk factors so that the new real scalar fields have canonical kinetic terms. At the end, finding the 
moduli mass squared eigenvalues boils down to diagonalizing the following matrix 

Mil = l^^^Ul^V^nUnl . (161) 

Unlike most of the other masses, the detailed form of the moduli mass matrix does depend upon the 
detailed form of Vx- Therefore we have resorted to numerical analyses in this case and found that there 
is one heavy modulus whose mass mainly depends on Q and for Q = 30 

M ~ 0(200 - 300) X m3/2 , (162) 

and lighter moduli with masses 



mj ~ 0(1) X m3/2 , i = l,N. (163) 

The heavy modulus arises from the fluctuation which deforms the volume of the three-cycle Vq, while 
light moduli originate from the fluctuations approximately preserving the volume and tangential to 
the hyperplane defined by 

N-s-Vq = 0. (164) 



The remaining light modulus represents the fluctuations of the hidden sector meson mixed with the 
geometric moduli. 
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C. Gaugino masses 



The universal tree-level contribution to the gaugino masses can be computed from the standard 
supergravity formula |24i] 



tree _ (-]Rr.\ 



where the visible sector gauge kinetic function is another integer combination of the moduli 



N 



i=l 



Note that, since the dominant F-term is that of the meson field, the gaugino masses at tree level will 
be suppressed wrt the gravitino mass. Since the scalar masses typically get contributions of order 
the expectation is to have light gauginos and heavier scalars, as we will indeed verify shortly. 

Plugging the solution for a ()124p into ()136|) while using the definitions (jllip we obtain 



^^/^i^'--#^(l + 77T^^ 1-3/2, (167) 



2si f 2Vx 

where we dropped the overall phase factor e"*'''^. It is now straightforward to compute the tree-level 
gaugino mass 

tree . . ^ f . . ^.Vx 



It is interesting to note that this formula is identical to the leading order expression previously obtained 
in [l| when one replaces the combination (j)^/Vx by the canonically normalized meson field. Here, again 
the suppression coefficient is completely independent of the number of moduli N as well as the integers Ni 
(AT?"*) appearing inside either the hidden sector ([76]) or the visible sector (I166P gauge kinetic functions. 
Moreover, all the detailed dependence on the individual moduli is completely buried inside the volume 
Vx and the gravitino mass 7713/2 (which also depends on Vx) and therefore expression ()168p is universally 
valid for any G2 manifold that yields positive solutions of the system of equations in (j96p . Hence, despite 
the presence of a huge number of unknown microscopic parameters, the tree-level gaugino masses in (jl68p 
depend on very few of them. Moreover, when the cosmological constant is tuned to a small value and 
Q — P = 3, the gaugino mass suppression coefficient becomes completely fixed! Indeed, using (jl46p and 
(fWl) for Q - P = 3 we obtain 

^/T ~ -0.0307 X m3/2 . (169) 

This result gets slightly corrected by the threshold corrections to the gauge kinetic function from the 
Kaluza-Klein modes computed in [12] 

1 ^ 



In the above formula, 71, is a topological invariant (Ray-Singer torsion) 



% = \n{Asm'^{57rw/q)) , 



(171) 
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where w and q are integers such that 5w is not divisible by q. In this case, the tree-level gaugino mass is 
given by 



m 



T/2 ~ -0.0307 r? X m3/2: 



where 



1 _ ^ 9gut q- 
87r2 



(172) 



(173) 



D. Anomaly mediated contribution to the gaugino masses 

Because of the substantial suppression of the universal tree-level gaugino mass, it makes sense to 
take into account the anomaly mediated contributions which appear at one-loop. The anomaly mediated 
contributions are given by the following general expression 25l | 



AM 



167r2 



\ a / \ a / a 

(174 

where Ca and are the quadratic Casimirs of the a-th gauge group and Ka are eigenvalues of 

the Kahler metric for the visible sector fields (j72p . Assuming the MSSM particle content, we have the 
following values for the Casimirs 
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c/(i): c„ = o Y.^: = -^ 

a 

SU{2): Ca = 2 ^C: = 7 
a 

5C/(3): Ca = 3 ^C° = 6. 



(175) 



Plugging the solution for a ()124p into (jl37p while using the definitions (|llip we obtain 

where we dropped the overall phase factor e'"^^^ . Combining (jl67p , (jl76p and using ([78|) , ([72]) , ([73]) and 
(|6ip we now compute the contributions 

.-."a,.>„;f„4(:.^) (e,.,|.Ij^)™3.. 

In the above we also used (|83p and (|84p together with the definition of Oj in ([16p as well as its contraction 
property (|17p . We also dropped unknown subleading contributions proportional to e^^'^F'^dic{si) ~ 
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Using the definition (jl74p we then obtain the following expression for the anomaly mediated contri- 
butions to the gaugino masses 



m 



AM 




(178) 



where we have explicitly separated the conformal anomaly contribution from the Konishi anomaly term 
using (1731) . 

Notice the appearance of the function c{si) which controls the size of the higher order corrections 
to the matter Kahler potential. As expected, when A = the Konishi anomaly vanishes in the exactly 
sequestered case [26|], i.e. when c{si) = 1. Again, in the leading order in l/Pg//, when c(sj) = the 
result obtained above is almost the same as the one in [l|. Just like in the case with tree- level gaugino 
masses, the above result is completely independent of the detailed moduli dependence of the volume Vx 
and therefore is completely general. In what follows, we will regard the value of the function c(sj) at the 
minimum of the scalar potential as a phenomenological parameter 



C[Si 



(179) 



When we set A = and Q — P = 3, tune the leading contribution to the vacuum energy by imposing 
the constraint ()146p . use ()147p and combine the above formula with the tree- level contribution ()172p . we 
obtain the following expression for the total gaugino masses 



-0.0307?? acuT 0.0364 3Ca 



a / 



+ 0.0749 (l-c)^Q 



(180) 



Note that as was previously pointed out in [26j], in the limit when c — )• 1 we obtain a particular type 
of a mirage pattern for gaugino masses 27|. However, as we will see below, in this limit the scalars 
become tachyonic and therefore, the exact mirage pattern is disfavored. An exact numerical computation 
confirms the above result giving 



Ma 



-0.03156 7? -Facc/T ^0.034086 ^3(7^-^(7°^ 0.07926 (1 - c) J]] C„"^ 



X m 



■3/2 



:i8i) 



Substituting the MSSM Casimirs p75]) into we then obtain 

Ml (-0.03156 r/ + ugut (-0.22497 + 0.52313 (1 - c))) x mg/a 
Ms K, (-0.03156 r/ + acuT (-0.03409 + 0.55483 (1 - c))) x mg/s 
M3 (-0.03156 r? + acuT (0.10226 + 0.47557 (1 - c))) x 7713/2 . 



(182) 



The form of (I182p allows us to see explicitly that for c = the Konishi anomaly contribution is larger than 
the contribution from the conformal anomaly by a factor of a few, which is what made the gaugino mass 
spectrum in Q very different from other known patterns. However, as we will see below, suppressing the 
scalar masses relative to the gravitino mass by tuning the coefficient c will automatically result in a large 
suppression of the Konishi anomaly. 
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E. Scalar s 



The masses of the unnormahzed scalars can be computed from the fohowing general expression 2j] 



^'Ib = H/2 + ^o) K^-p - e^F^F'^ [drndnK^p - d^K^^K^^dnKs^) . (183) 

Since the vacuum energy is tuned to zero we set Vq = in the above. Using ([68]) . ([78]) . ([83]). ([8^ . the 
contraction properties ([8]), ([9]), ([6T]) and the F-terms ()167p . ()176p we obtain from ()183p in the leading 
order 

mf^ « (1 - c(.,)) - ^ (mf/f ) K,^ + A (mf/f ) ' K,^ , (184) 

where for consistency reasons we only kept contributions linear in c{si) and dropped unknown subleading 

terms proportional to the derivatives of c(sj), such as e.g. F^F^ didjc{si) ~ (^^1/2^ SiSjdidjc{si) . In 
the above derivation we also used the following properties 

e^rPd^ch^k = 7 (mf/-)' ^ e'^F-F^dndrn^ = mlj^ - ^ (mf/;)' . (185) 

Notice that despite the presence of the derivatives of the Kahler metric K^^ in the definition ()183p , the 
final expression (11841) contains K^j^ only as an overall multiplicative factor. This happened because the 
moduli F-terms (11671) up to a phase are essentially given by = 2si X m}y2 s-nd the matrix Kapi^i) 

is a homogeneous function satisfying (I6ip . Therefore, diagonalization and canonical normalization of 
the corresponding kinetic terms automatically results in universal masses for the canonically normalized 
scalars 



mt « 1 



c(..)) (m^/, - I (<-) + A (mf/-) ' . (186) 



After setting A = 0, we tune the leading contribution to the vacuum energy by imposing the constraint 
(fTi5|) and use (fTM|) to obtain from (fTM|) 

m„ w (1 - 0)^/^0.999 m3/2 , (187) 

where we again treat the value of the function c(sj) for a given vacuum as a phenomenological parameter 
(|179p . A numerical computation in this case gives excellent agreement 

ma -{I- c)i/2o.998 ^3/2 w (1 - c)i/2 ^^^2 ^ (188) 

Again, for c = we recover the old result in [l| where all the scalars have a flavor-universal mass equal 
to the gravitino mass. 

Furthermore, the anomaly contributions to the scalar mass squareds are suppressed relative to the 
gravitino mass and since we wish to consider generic 0{1) values of (1 — c) we will neglect such contri- 
butions. Concretely we are going to consider only those values of < c < 1 which give 

1 m„ 
levr^ m3/2 

such that the anomaly mediated contributions to the scalar masses can be safely neglected. However, 
one can certainly extend our model and include such contributions. Once again, the result above is 
completely independent of the details of Vx and therefore holds for any G2 manifold that solves the 
system (I96p with Oj > such that the Kahler metric at the minimum is positive definite. 
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F. Trilinear couplings 



The unnormalized trilinear couplings for the visible sector fields can be computed from the following 
general expression [23] 

= J^e"^^'^" + ^raYL^ - (k^'Pdmkp^YL + (a O /3) + (o O 7)) 1 , (190) 



\w\ 



where {a , /3 , 7} label visible sector matter fields and Y^^^ are the unnormalized Yukawas that appear 
in the superpotential. Recall that the Yukawa couplings 1^^^ arise from the membrane instantons wrap- 
ping associative cycles Q'^^'^ , which connect isolated singularities supporting the corresponding matter 
multiplets. They are given by 

Y^^^ = Co.p^e'^^^^<'"'^. (191) 

The integer combination of the moduli Vqafj-, = m"^'^ Si gives the volume of the associative cycle Q"'^''' 
connecting co-dimension seven singularities a, (3 and 7 where the chiral multiplets are localized. The 
coefficients Cq,^^ are constants. The relation between the physical and unnormalized Yukawa couplings 
is given by 

Y^p^ = ^e^/2y^^^ {K^kpk,) '^'^ . (192) 
Using (I167p and ()19ip we can compute the contribution 

e^l'F-dmY!^,^ = YLp,^^ (1 + (Q^^) yQ'^f^^^m ■ (193) 

Similarly, using §Ei), ([78]), §3\), ([HID, the contraction properties §1, (HD, §B) and the F-terms (fWl) . 
(flTHIl we find 

.*v.^™^*.a„A-^ 4 (1 . ^i,) (c,.)| . ^) ... . im 

where for consistency reasons we only retained contributions linear in c(sj) and dropped unknown sub- 
leading terms proportional to e^/"^ F^dic{si) ~ m^y^Sidic{si). Again, in the above expressions we did not 

display the overall phase factor e~^'^^ . Using the definition (I190p along with (I177p . (I193P and (jl94p we 
obtain the following expression for the physical (normalized) trilinear couplings at tree-level 



A.,-, = ( 1 + (d - + ] ™3/. , (1%) 



which gets reduced to the result in [1] when c = and A = 0. Once again, the detailed structure of 
the volume Vx played absolutely no role in our ability to obtain the above expression for the tree-level 
trilinear couplings. The actual volumes of three cycles Vqap-y do depend on the microscopic properties 
of G2 manifolds and in our general framework these parameters remain undetermined. However, below 
we will present a good argument for dropping such volume contributions completely when the third 
generation trilinear couplings are computed. 

Setting A = and Q — P = 3, when the leading contribution to the vacuum energy is tuned we obtain 
for the reduced trilinears (the physical trilinears divided by the physical Yukawa couplings) 



A^p^ = ^ = (1.41(1 - c) + 0.386 X Vq.,-,) ^3/2 . (196) 
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Prom the corresponding numerical calculation we obtain the following result 

i„/37 = (1.494(1 - c) + 0.3966 x Vgcfi-y) ■ (197) 

Since the physical Yukawa couplings for the third generation fermions are much larger than the first 
two generation Yukawas, one can typically neglect the trilinears for the first and second generations. 
Moreover, the large size of the third generation Yukawas implies that the volumes of the three-cycles 
of the corresponding membrane instantons are very small. In fact, because the top Yukawa is of order 
one, one can assume that the point pi supporting the up-type Higgs 5 of SU{5) coincides with the point 
Do supporting the third generation 10, so that the coupling ff^lOslOs has no exponential suppression 
[5(], [12l |. At the same time the point supporting the down- type 5 Higgs and the point p4 supporting 
the third generation matter 5 are distinct but still close to P2,Pi so that the coupling of HdbslOs 
which accounts for the bottom(tau) Yukawa is slightly smaller than the top Yukawa at the GUT scale. 
These considerations completely justify dropping the corresponding Vqap-/ terms for the third generation 
trilinears which then become simplified 

At = Ab = Ar^ 1.494(1 - c) my2 ■ (198) 

For generic values of c the trilinears are of the same order as the gravitino mass. In the limit c — )• 1, 
the reduced trilinear couplings at tree-level become suppressed relative to the gravitino mass. Note that 
as c approaches one, the suppression of the trilinear couplings above is much stronger than that of the 
scalars. In this case, the anomaly-mediated contributions may become com par able to the tree-level ones 



and therefore must be taken into account. General expressions given in 28|] can be simplified in the 
nearly sequestered limit as 

Ar = - Y^7. [e-^'W - yi-F-K:^ + ii^X. m,i2 , (199) 

where the last term denotes the unknown contributions vanishing in the sequestered limit. Note that 
such terms are suppressed compared to the tree- level piece (|198p due to the loop factor. As long as (1 — c) 
is small enough, they become subleading and we will drop them in further analysis. Using (jl77p and 
substituting the corresponding MSSM expressions for 7^8, where we set gi = 92 = 93 = Qgut, we obtain 
the following expressions for the anomaly mediated contributions to the reduced trilinear couplings 

^t" » {-U.r . + f 1 - f # + ^) ) "3. (.00) 



16^2 ^ "^t -V V 3 V {Q-P)(t>lJ \Vx Peff, 



-\AM 



When we set Q — P = 3, tune the tree-level vacuum energy by imposing the constraint (I146p . use (1147(1 
and combine the above formula with the tree- level contribution ()197p . we obtain 

At « 1.41(1 - c)m3/2 - 0.0029 {-^qIut + + n') "^3/2 (201) 

A, « 1.41(1 - c)m3/2 - 0.0029 (-y^ci/T + Yt + 6^' + "^3/2 

Ar « 1.41(1 - c)m3/2 - 0.0029 (-y^GC/T + + 41;') "^3/2 • 
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Numerical computations give the following expressions for the total reduced trilinears 

46 

44 





i 1.494(1 


- c)m3/2 




a 1.494(1 


- c)"i3/2 




5i 1.494(1 


- c)m3/2 



(202) 



24 o n 

■y5GC/T + 3yb' + 4y 

which demonstrate a fairly high accuracy of the analytically derived result in (120ip . 

G. /i and -terms 

The full hidden sector plus visible sector Kahler potential and superpotential can be written in the 
following general form 



i^total = i^(s*,</',<^) +^«^(Si,</>,<^)Q"gt'' + Z„/3(Si,</>,<^)Q"Q''+ /i.e. (203) 

4p + Ai'Q^Q'' + Y!,p^i 



W = Wnp + fi'Q'^Q" + Yl^g^D'^Q^Q^ + .... 



Here, (j) denote the hidden sector matter fields while are visible sector chiral matter fields where 
K^^{si,(j),(j)) is the visible sector Kahler metric and Y^j^^ are the corresponding unnormalized Yukawa 
couplings. It can be shown that the supersymmetric mass parameter can be forbidden by requiring 
certain discrete symmetries which are also used in order to solve the problem of doublet-triplet splitting 
[igl - Hence, in our analysis we will rely on the Giudice-Masiero mechanism [i^] in generating effective /i 
and Bfi terms where the bilinear coefficient Z^{si, (p, (p) in (j203p plays a key role. The general expressions 
for the normalized /x and Bfi are given by [2J] 



(^e^/V + my,Z - e^/^F-a^z) (K^^K^J"^/^ (204) 



\Wnp 



+ (2m2/2 + Vo)Z- m3/2e^/2^'^a^Z + m3/2e^/2^"(a„Z - HKh^Kh,)) 



-1/2 



where we can set = 0. Unfortunately, at this point we do not have a reliable way to compute the Higgs 
bilinear Za/3{si, (p, 4>) for G2 compactifications. Therefore, in our analysis we will parameterize the /x and 
B/i terms as follows 

fi = Z^ffmy2 (205) 
Bfi = Zgjj 7713^2 ! 

and treat Z^^^ and Z^^j as phenomenological parameters. Naturally, we expect that Z^^j ~ C'(l) and, 
as we will see in the next section, tuning n parameter in order to get the correct value of the Z - boson 

1 2 

mass boils down to tuning the values of Z^jj: . 
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VI. GENERALIZATION TO THE CASE WHEN k{s,) IS A NON-TRIVIAL HOMOGENEOUS 

FUNCTION 



Recall that the above results have been obtained assuming that the factor appearing in the 

Kahler potential (|78p of the hidden sector matter is a pure constant. In this section we briefly outline 
the main results for the case when K{si) is a general homogeneous function of the moduli of degree zero, 
satisfying the property (|3ip . Here we will not give any explicit analytic derivations (these are quite tedious 
and would mostly resemble the computations in the preceding sections) and instead present numerical 
evidence that most of results obtained for the simplified case K{si) = 1 can be directly extended to the 
more general scenario. 

The main difference from the previous case is in the form of the moduli vevs at the minimum of the 
scalar potential. These are now given by 

iV, r * 3 V UniQ-P) ' 

where r « 3/2 when Q — P = 3. Note that the vev (/>^ of the canonically normalized effective meson field 
at the minimum is given by the same expression as in the case when K{si) = 1: 

where Peff is exactly the same as in (I14ip . Keep in mind that the analytic expression (I207p is only valid 
when the leading contribution to the vacuum energy is tuned to zero. Parameters Cj are defined as 

din k( Si) 

a = Si , no sum over i , (208) 

dsi 

and satisfy 

TV 



^Ci = 0, (209) 



=1 



because K{si) is a homogeneous function of degree zero. At the minimum of the potential parameters 
and Ci can be determined by solving a system of 2N coupled transcendental equations 



dK 

'' dsi 



+ 3ai =0, (210) 



5 In nisi) 

s. 



dsi 



Once again, the structure of the soft supersymmetry breaking terms remains virtually unchanged in 
the leading order in l/Pg// expansion and does not depend on the precise details of the function K{si). 
In fact, the only modification compared to the previously derived expressions is the replacement of the 
following combination 
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whose vev at the minumum is given by (|207p and is exactly the same as before! However, it turns out 
that the soft super symmetry breaking terms now become shghtly sensitive to the compactification details 
via the subleading corrections. The most sensitive parameter is the tree- level gaugino mass. Up to an 
overall phase it is given by 

tree^ "^3/2/, , 2Vx , M ^ _ "^3/2 A , 2 5 



"^''^ " Peff V^{Q- P)-{S^)K + ^ \Peff)) " Peff V ^ {Q - PW. ^ Peff) ' ^'''^ 

where we introduced a phenomenological quantity 5 ~ in order to parametrize the additional 

correction. In the numerical toy examples we studied, we obtain a 0(1 — 10)% variation in the value of 
the tree-level gaugino mass, while the other soft terms vary by less than 1%. For the sake of completeness 
we shall list the expressions for the remaining soft breaking terms 



AM _ QGt/T 



(^3Ca - c:^ ^1 + x; c:^ x m,/^ , (213) 



2 ) 



« fl 



K,^l-Ul+ ^^ + (214) 



where we defined 



3 V (Q-P)K(si)0gy V VX Peff, 

To illustrate the high accuracy of the analytical results presented here we present a simple toy example 
with two moduli. 

2 

yx = sxsT. k(s*) = 1 + 4- (215) 

•^2 

For the following choice of the parameters in the superpotential 

^1 = 27 , ^2 ~ 2.638 , P = 27 , Q = 30 , iVi = 1 , A^2 = 2 , (216) 

where was tuned to cancel the leading contribution to the vacuum energy, we obtain numerically 

sif« 71.67, S2f« 13.02, 00^^7.24, ^ 0^5:^ 0.746, Pg// 61.68 . (217) 

Notice that the above values of and Pg// are extremely close to those in (jl47p and (|146p . obtained 
numerically for the case when K(sj) = 1. These values are also in very good agreement with the corre- 
sponding analytical results. By solving the system (j210p for the above example we obtain 

oi = 1 , 02 = - , ci = -C2 w 1.939 , (218) 
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and using the above values in the analytic expression ()206p with Pg// ^ 61.68 we find 

si ^ 72.5, S2 ~ 12.8, (219) 

which agree well with the numerically obtained values in (j217p . To verify the tree- level gaugino mass 
formula numerically we need to know the integers A'';^''^ and A^J'*^ for the visible sector gauge kinetic 
function. Here we list three representative examples where we varied N^^^ and N2^^, while keeping 
everything else fixed 

> m\j^ ^ -0.029 X m3/2 , (220) 

> mflf w -0.032 X m3/2 , 

> m\j2^ « -0.037 X m3/2 , 

which demonstrate a mild dependence of the tree-level gaugino mass on the compactification-specific 
details. Using 0^ w 0.746 together with Peff ~ 61.68 in the analytic formula (j212p for Q — P = 3 we 
obtain 





= 1, 


Nr 


= 




= 1, 




= 3 




= 0, 


Nt 


= 1 



m 



tree 
1/2 



(0.031 + 0.00026 x 5)x mg/g , (221) 



which is in fairly good agreement with the numerical results. The numerical results for the remaining 
soft terms are given by the following expressions 



m^^^ w acuT ^0.0337 (^3Ca - J2 + O-O'^^^ (1 - c) C^^ x my2 , (222) 
Aa « 1.49(1 - c)m3/2 - 7a X 0.0027 x mg/a 
ma w (1 - c)^/^ m3/2 , 

Kiss 0.424, /Ca ~ 1-494(1 -c) , 

and are virtually unchanged compared to the numerical results in (jl8ip . (j202p and (|188p . computed for 
the case when n^Si) = 1. The values obtained from the corresponding analytic expressions are in good 
agreement with the numerical values above and give essentially the same results as in the case when 
Ac(sj) = 1 because the values of Peff and cf^ barely changed. Thus, the effect of including a non-trivial 
function K(sj) in the Kahler potential for the hidden sector matter fields can be reliably described by a 
single parameter 6 that appears in the subleading contributions to the tree-level gaugino mass, while the 
remaining soft terms stay essentially unaffected. 



VII. ELECTROWEAK SCALE SPECTRUM 

In order to obtain the corresponding MSSM spectrum at the electroweak scale we need to RG-evolve 
all the masses and couplings from the GUT scale down to the electroweak scale. This procedure was 
described in great detail in Here we will only highlight a few important points and give the final 
results. 

As we have seen in the previous section, at the GUT scale, the gaugino masses are non-universal and 
highly suppressed relative to the gravitino mass. On the other hand, unless c is very close to one, the 
scalars, trilinear couplings and the /i-term are all of order Hence, we can define a scale at 
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which ah the heavy states decouple and the effective theory below that scale is the Standard Model plus 
gauginos. More specifically, we can choose the decoupling scale nis to be the geometric mean of the stop 
masses 

nis = ^mimi^ . (223) 

This is okay as long as the mass differences between the lightest stop and the other heavy states is not too 
large. Then, the running can be done at one loop in two stages with tree- level matching at the scale m^. 
This method, however, does not capture the two-loop effects, which may give significant contributions to 
the running. Thus, in what follows we will utilize the SOFTSUSY package ^] and perform the running 
at two-loops with the full the MSSM spectrum and accound for the effects from the heavy scalars via 
threshold corrections. 



A. Gauginos 

As one notices from (j222p . due to the anomaly mediated contribution, the gaugino masses are sensitive 
to the value of ugut- However, the value of acuT is only determined once we know the exact spectrum 
and run the gauge coupling up to the GUT scale. Therefore, there is a feedback mechanism, which allows 
us to completely fix the gaugino masses by imposing the gauge coupling unification. In practice, we first 
pick an initial value of acuT ~ 1/25, compute the gaugino masses, scalar masses, trilinears, etc. at the 
GUT scale and run them down to the electroweak scale where we compute the spectrum. We then run 
the gauge couplings up using two-loop RGEs to check if they unify at the same value of acuT as we 
chose to compute the gaugino masses. If there is disagreement, we change the value of acuT by a small 
increment and repeat the steps until there is a match. In addition, parameter r] which appears inside 
the gaugino masses and was defined in (jl73p can be safely set to one. This is because as one varies the 
integers w and q inside (117ip over a reasoble range, the torsion, unless specifically tuned, is so small that 
that the KK threshold corrections can be neglected. 

Since Mniggsino ~ ~ 0(7713/2), there is a substantial threshold contribution from the Higgs- Higg sino 



loops which has to be taken into account when computing bino (Mi) and wino (M2) masses [4], 32l434i|: 
\T\i ~ «i,2 /isin(2/3) /i^ ai,2 ai,2^i 

AMi 2 ~ — T^-/ ^ In — T ^ — ; — M = — H-^p/-f"^3/2 • (224) 



,2 



In the above expression we expanded the logarithm using ~ 1 and used tan/3 ~ C(l)- The latter is 

especially true when 2 Z^^j ^eff relied on the fact that the supersymmetric /u-term almost 

does not change with the RG evolution so one can use (|205p . Since 777.3/2 ~ C'(10)TeV, the above 
correction to M2 can be as large as a few hunrded GeV. It turns out that this contribution is intimately 
related to the value of parameter c that directly affects the scalar masses and indirectly forces the value 
of /X to get smaller, as the scalars get lighter. For < c ^ 0.05 and 0.8 ^ /x/777,3/2 the contribution (I224p 
is actually large enough to completely alter the nature of the LSP depending on the sign of ^. In this 
respect, the sign of 5 that parametrizes the sub leading corrections to the tree- level gaugino mass ()22ip 
also plays an important role. In particular, from the left plot in Fig [1] where we picked c = there is a 
region where 5 ^ —0.7 such that the LSP is Wino-like, while for 5 ^ —0.7 the LSP becomes mostly Bino. 
Furthermore, as one can see from the plot, there exists a small range of values where Mi and M2 become 
nearly degenerate. This is certainly an intriguing possibility, which may provide for a well-tempered 
neutralino candidate [s^]. Note that in the Wino-like LSP case, the lightest chargino and neutralino are 
degenerate at tree-level, i.e. Xi = Xi = M2. However as we take into account the 1-loop contribution 



from the gauge bosons 3j], this degeneracy is removed, as is seen from the corresponding entries in Table 
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II. Such splitting was discussed in detail for the pure anomaly mediation scenario in [35|, |36| and is given 

by 



^M.^ioo, = ^[f[j^)--wf\j^^]- swf (0) 



(225) 



where /(a) = Jq dx{2 + 2x)ln[x'^ + (1 + x)a'^]. Typically we obtain 160 MeV < AMi^ioop < 200 MeV. 
Because of this, the lightest charginos are quazistable and decay into LSP plus soft pions or soft leptons. 
In the collider context such decays would take place well inside the detector leaving short charged tracks. 
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FIG. 1. Left panel: Gaugino mass parameters at the electroweak scale as functions of 5. In the above computation, we 
used SOFTSUSY with the high scale input 7713/2 ~ 30TeV, /i < 0, tan/3 = 2.6, c = 0. We have verified by examining the 
neutralino mixing matrix that at 5 ~ —0.7 the LSP type changes from Wino to Bino as S is increased. 

Right panel: Gaugino mass parameters at the electroweak scale as functions of 7713/2- For each data point we checked the 
neutralino mixing matrix to confirm that the LSP is Wino-like. The plot was generated using SOFTSUSY with /i < 0, 
c = 0, S — —5, tan /3 = 2.6. The KK threshold correction to the visible sector gauge kinetic function was neglected in both 
plots. 



In addition to (|224[) . the EW threshold corrections from gaugino-gauge-boson loops must also be 
included, especially for the gluino 

AMf^i = ^ (3 In + 5) M3 . (226) 

Unfortunately, we have no technical handle on the size of the parameter c, though we expect it to 
be small and hence a Wino LSP is quite generic when c ~ 0. As we increase the value of parameter 
c, the LSP quickly becomes Bino-like. There are two reasons for this effect. First, the ratio M2/M1 at 
the GUT scale grows as c is increased from zero to one. At the same time, the scalars and higgsinos 
become lighter relative to the gravitino mass. In particular, for a fixed the lower bound on the 

Higgs mass forces us to consider somewhat larger values of tan (3 which in turn leads to smaller values 
of fj, thus significantly reducing the contributions (12240 from higgsinos. Recall that it was primarily due 
to this contribution from heavy higgsinos for fj, < that the LSP could become Wino-like. Thus, due to 
the increase of M2/M1 at the GUT scale and the decrease in the higgsino mass, the Wino-like LSP case 
becomes rapidly excluded as we increase the value of c. Benchmarks 4 and 5 in Table II demonstrate 
that for generic values of c > the LSP is always Bino-like. 



Of course, pure Bino LSP is almost certainly excluded by the standard cosmological considerations j33l |. 
Namely, because binos do not annihilate efficiently, the dark matter relic density becomes unacceptably 
large. However, this problem can be avoided when the higgsinos, which annihilate efficiently, are light 
enough to mix with gauginos. If the higgsino component of the LSP is significant, it can easily reduce 
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the relic density to acceptable levels by increasing the annihilation crossection of the LSPs. It turns out 
that for generic values of c, < c < 1, the higgsinos are always much heavier than the gauginos. This 
is because at the decoupling scale m^, the yU^-term must be of the same order of magnitude as mjj to 
give a correct value of the Z-boson mass, and since for typical values of c we get >> Mf 2: the 

higgsinos do not mix with gauginos. 

The only way to make the higgsinos light is to consider values of c very close to one. In this limit 
both scalars and higgsinos also become highly suppressed relative to the gravitino mass. As we will see 
in the analysis of this limit, the higgsinos can become as light as the gauginos and for certain values of 
tan /3 can give a well-tempered LSP candidate. 



B. Squarks and sleptons 

Recall that at the GUT scale all the squarks and sleptons have a universal mass ()188p . which for 
generic values of c (0 < c < 1) is smaller but nevertheless typically of the same order of magnitude as the 
gravitino mass. However, as we evolve these down to the electroweak scale, the third generation scalars 
become significantly lighter whereas the first and second generation scalars experience a very mild change 
in their masses. Indeed, because the third generation Yukawa couplings are large, the stops, sbottoms, 
and staus are affected through the corresponding trilinear couplings (1202 p . which are of 0(771,3/2). As 
one can see from Table II, this effect is especially dramatic for the lightest stop ti. Yet, it is still much 
heavier than the gauginos and is effectively decoupled from the spectrum at the electroweak scale. 

However, since gluinos can be pair produced at the LHC via gluon fusion, the gluinos (which have to 
decay via a quark-squark pair) have a sizeable branching fraction into top-stop - precisely because the 
stop is the lightest squark. This leads to events containing up to four top quarks at the LHC [s^ 



C. Radiative electroweak symmetry breaking 

The existence of the electroweak symmetry breaking (EWSB) in the effective theory below the decou- 
pling scale rus is determined by whether there exists a negative eigenvalue in the Higgs mass matrix 



(227) 



at the scale where the scalars decouple [37|]. Recall that at the GUT scale, rn^H^ = t^'h^ = (1 ~ c) rn^i^ 
whereas ^ = Z^j^m^i2 and = Z^^^m^j^- It is well known that the positive contribution into the 
running of the up Higgs mass parameter squared from the stop is crucial for radiative EWSB as it drives 
Tn\^ negative. It turns out that for a fixed value of the gravitino mass 7773/2, as we vary parameter c 
there exists a narrow range of values ^e// which the matrix ()227p has a negative eigenvalue above 
the decoupling scale 777^ defined in (j223p . 

However, unless we force -^e// << 2 Z^^^, all the entries in the above matrix are 0(7773/2)- Therefore, 
both the lightest Higgs mass and the Z-boson mass naturally come out to be of 0(7773/2). For the same 
reason, there is little mixing between the two Higgs doublets and naturally tan/3 ~ 0(1) is predicted 

1 2 

in our framework. In practice, parameters Z^'j^ must be tuned in such a way that the corresponding 



eigenvalue turns negative right at the decoupling scale [37| so that 777^ ~ 91 GeV. This fine tuning is 
a manifestation of the so-called little hierarchy problem - the hierarchy between the electroweak scale 
Mew ~ 0(100) GeV and the scale where the scalars decouple 777^ ~ 0(10) TeV. Once Mz is tuned, the 
Standard Model Higgs mass turns out to be rrih < 130 GeV. 
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TABLE II. Low scale spectra for seven benchmark G2-MSSM models generated by SOFTSUSY package. All masses 
are in GeV. The top mass was taken to be mt = 171.3 GeV. Here we only give absolute values of the gaugino 
masses and suppress the relative phases. The; spectra are largely determined by four parameters rnz/2, S, c and 
tan/3. The Kaluza-Klein threshold corrections to the gaugino masses have been neglected. For the above spectra, 
the gauge couplings unify at the value of aQ^j, « 24.2 at the scale Mqut » 2 x 10-^^ GeV. 
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D. Approximately sequestered limit 

For completeness, in this subsection we are going to examine the spectrum in the Umit when parameter 
c approaches one, though this is presumably not a case which really originates from a "typical G2- 
manifold". Since the scalars will now be light, the full MSSM running must be implemented from the 
GUT scale to the electroweak scale. To generate the sparticle spectrum at the electroweak scale we utilize 
the SOFTSUSY package [sij with the GUT scale soft terms specific to our model. Before we move to 
computing the MSSM spectrum at the electroweak scale, we would like to point out a useful relation 
between the soft terms we have computed. Namely, we can use the expression for the scalar masses (I188p 
to express parameter c as 

c=l-(j^\\ (228) 

We can then use (I228P to recast the gaugino masses and reduced trilinears as follows 

2 

Ml ^ I -(0.031 + 0.00026 x 5) ?? + acuT ( -0.22497 + 0.52313 



m. 



M2 ^ |^-(0.031 + 0.00026 X 6)7] + acuT |^-0.03409 + 0.55483 
M3 ^ \ -(0.031 + 0.00026 x 6)i] + acuT ( 0.10226 + 0.47557 



m, 




m, 



m 



At « 1.494 (^-^^ m3/2 - 0.0027 (-y^GC/T + 61"*' + n') x m^/^ (230) 
Ab « 1.494 (^) 7713/2 - 0.0027 (-^ghuT + + ^^h + >f ) x ms/a 



3/2 

A, « 1.494 I m3/2 - 0.0027 {-^-qIut + 3^' + 4lf ) x ^3/2 . 



"I3/2 

From the above expressions it is clear that as the value of c approaches one, the suppression of the scalars 
at the GUT scale relative to the gravitino mass results in a much larger suppression of both the Konishi 
anomaly contribution inside the gaugino masses and the tree- level term in the reduced trilinear couplings. 
Such a non-trivial interplay among the soft terms is very remarkable and leads to some very interesting 
MSSM spectra unique to our model. 

For convenience, instead of parameter c we choose to use the universal scalar mass rn-Q, as one of the 
inputs. Because we cannot yet compute the /x and B[i terms from first principles, we use tan /3 as an 
input and tune [i to get the correct value of the Z boson mass. Although acuT {Qgut) appears as a 
parameter, it is actually fixed at the value where the gauge couplings unify acuT ~ 1/25. The GUT 
scale physical Yukawa couplings in (j230p are fixed once tan (3 is chosen. In total, we have five relevant 
high scale input parameters: 

"^3/2 ; n^a, 1], 5, tan/3. (231) 

However, as was mentioned earlier, parameter r] which appears inside the gaugino masses can be safely 
set to one. As we vary the ratio ma/?n,3/2, we need to make sure that we are in the regime where the 
anomaly mediated contributions to the scalar masses can still be neglected. 
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3192.1 




914.9 


1194.5 


1193.9 


1477.8 


2527.2 


3225.6 




733.7 


992.0 


1032.1 


1193.2 


2060.8 


2813.1 




714.8 


1003.1 


1026.6 


1164.8 


2004.4 


2781.8 




809.2 


835.9 


836.0 


1223.4 


2048.4 


2134.5 




755.4 


816.3 


829.7 


1152.5 


1940.2 


2159.6 


mho 


111.4 


115.0 


115.1 


115.7 


119.3 


121.7 


mHo 


426.5 


444.6 


637.3 


653.3 


1140.7 


1774.5 


rriAo 


426.4 


444.6 


637.2 


653.2 


1140.6 


1774.5 


mH± 


435.4 


452.6 


642.4 


658.4 


1143.4 


1775.8 


At 


445.5 


810.5 


830.3 


829.5 


1485.9 


2397.4 


h 


561.2 


954.2 


1079.2 


1014.4 


1811.7 


3112.4 


K 


159.9 


192.7 


272.2 


257.2 


473.0 


864.9 



TABLE III. Low scale spectra for six benchmark G2-MSSM models with approximate sequestering generated by 
SOFTSUSY package. All masses are in GeV. The top mass was taken to be rrit — 172.3 GeV. The thermal 
relic density was computed by combining the output from SOFTSUSY with micrOMEGAs software. For all the 
benchmarks we fixed S = 2, rj — 1 and varied parameters nia, 1713/2 and tan/3. For the above spectra, the gauge 
couplings unify at the value of acuT ~ 1/25. 



46 



Table III contains six benchmark points generated by SOFTSUSY. To compute the thermal relic 
density we used SOFTSUSY generated output as input for the micrOMEGAs program [38]. Note that 
the first benchmark point is experimentally excluded by the lower bound on the Higgs mass. As one can 
see from the entries for Mi and M2, the LSP is mostly Bino. Thus, for a generic value of tan (3 the thermal 
component of the relic density turns out to be much larger than the observed value Qch'^ ~ 0.1099ib0.0062 
[4o| . However, for large values of tan (3 there is a window where the Higgsino component of the LSP can 
be large enough to reduce the relic density. Hence, for each set of ma and we can tune tan /3 so 
that the value of Qch'^ is compatible with experiment. For example, Points 2 and 3 with rua = 800 GeV 
and = 16TeV both satisfy the experimental bound on the relic density but have different values 

of tan/3, i.e tan/3 = 43 and tan/3 = 38. In this case, if tan/3 is within a window 38 ^ tan/3 ^ 43, the 
thermal component of the relic density is too low and one may have to appeal to non-thermal mechanisms 
0]. Note that the SU{5) GUT scale relation Yh ^Yj- automatically holds in our construction. However, 
for large tan f3 the values of all third generation Yukawa couplings become comparable. Namely, for the 
benchmarks in Table III, we obtain 

0.524 < < 0.547 , 0.28 < Yb ^ Yr < 0.39 . (232) 
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